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ABSTRACT 

In  late  2008,  the  quasi-persistent  neutron  star  X-ray  transient  and  eclipsing  binary 
EXO  0748-676  started  a  transition  from  outburst  to  quiescence,  after  it  had  been 
actively  accreting  for  more  than  24  years.  In  a  previous  work,  we  discussed  Chan¬ 
dra  and  Swift  observations  obtained  during  the  first  five  months  after  this  tran¬ 
sition.  Here,  we  report  on  further  X-ray  observations  of  EXO  0748-676,  extending 
the  quiescent  monitoring  to  1.6  years.  Chandra  and  XMM-Newton  data  reveal  qui¬ 
escent  X-ray  spectra  composed  of  a  soft,  thermal  component  that  is  well- fitted  by  a 
neutron  star  atmosphere  model.  An  additional  hard  powerlaw  tail  is  detected  that 
changes  non-monotonically  over  time,  contributing  between  4  and  20  percent  to  the 
total  unabsorbed  0.5-10  keV  flux.  The  combined  set  of  Chandra ,  XMM-Newton  and 
Swift  data  reveals  that  the  thermal  bolometric  luminosity  fades  from  ~  1  x  1034  to 
6  x  1033  (d/7.4  kpc)2  erg  s_1,  whereas  the  inferred  neutron  star  effective  temperature 
decreases  from  ~  124  to  109  eV.  We  interpret  the  observed  decay  as  cooling  of  the 
neutron  star  crust  and  show  that  that  the  quiescent  lightcurve  of  EXO  0748-676  is 
markedly  shallower  than  that  observed  for  three  other  neutron  star  X-ray  binaries 
that  underwent  prolonged  accretion  outbursts. 

Key  words:  accretion,  accretion  disks  -  binaries:  eclipsing  -  stars:  individual 
(EXO  0748-676)  -  stars:  neutron  -  X-rays:  binaries 


1  INTRODUCTION 

EXO  0748-676  is  an  intensively  studied  low-mass  X-ray  bi¬ 
nary  that  was  initially  discovered  with  the  European  X- 
ray  Observatory  SATellite  ( EXOSAT )  in  1985  February 
(jParmar  et  al.l  1198511 .  In  retrospect  the  source  already  ap¬ 
peared  active  in  EXOSAT  slew  survey  observations  several 
times  beginning  1984  July  feevnolds  et  al .111 9 9 III ) .  The  de¬ 
tection  of  type-I  X-ray  bursts  fe.g..  iGottwald  et  al.l  1 1 986li 
identify  the  compact  primary  as  a  neutron  star.  A  few  X-ray 
bursts  have  been  observed  that  exhibited  pliotospheric  ra¬ 
dius  expansion,  which  allows  for  a  distance  estimate  towards 
the  source  of  D  =  5.0  —  8.3  kpc,  depending  on  the  compo¬ 
sition  of  the  photosphere  ll Wolff  et  al.ll2005:  Galloway  et  af.l 
2008).  The  system  exhibits  irregular  X-ray  dips  and  displays 

*  e-mail:  degenaar<Q!uva.iil 


eclipses  that  last  for  ~  8.3  min  and  recur  every  3.82  hr, 
which  allow  the  unambiguous  determination  of  the  orbital 
period  of  the  binary  llParmar  et  al.lfToSlil:  IWolff  et  ~al]l2009l') . 

At  the  time  of  its  discovery,  EXO  0748-676  was 
detected  at  2-10  keV  luminosities  of  ~  (1  —  7)  x 

1036  (D/7.4  kpc)2  erg  s_1  ((Parmar  et  al.l  1 1 9861).  How¬ 
ever,  several  years  prior  to  the  EXOSAT  detections,  it 
was  serendipitously  observed  with  the  EINSTEIN  satel¬ 
lite  in  1980  May,  displaying  a  0.5-10  keV  luminosity 
of  ~  5  x  1033  (D/7.4  kpc)2  erg  s_1  llParmar  et  al.l  1 19861: 
iGarcia  fe  Callananl  1 1 9 9 ill .  The  source  can  thus  be  clas¬ 
sified  as  a  transient  X-ray  binary.  Nevertheless,  such 
systems  typically  exhibit  accretion  outbursts  that  last 
only  weeks  to  months  (e.g.,  IChen  et  al.l  Il997ll.  whereas 
EXO  0748-676  was  persistently  detected  at  luminosities  of 
~  1036-37  (D/7.4  kpc)2  erg  s-1  by  various  satellites  for  over 
24  years.  Similar  prolonged  accretion  episodes  continuing 
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for  years  to  decades  have  been  observed  for  a  few  other  sys¬ 
tems,  which  are  termed  quasi-persistent  X-ray  binaries  (e.g., 
IWiinandslEooii. 

In  2008  August-September,  observations  with  the  Pro¬ 
portional  Counter  Array  (PCA)  onboard  the  Rossi  X- 
ray  Timing  Explorer  ( RXTE )  and  Swift’s  X-ray  Telescope 
(XRT)  indicated  that  the  X-ray  flux  of  EXO  0748-676  was 
declining  (IWolff  et  al.ll2008alfbh .  Optical  and  near-IR  obser¬ 
vations  of  the  optical  counterpart,  UY  Vol,  performed  in 
2008  October  showed  that  the  optical  emission  had  also 
faded  compared  to  the  brighter  X-ray  state  ()Hvnes  fe  Jonesl 
l200Sl.  [2003:  iTorres  et  al]l200Stl.  These  events  indicated  that 
the  accretion  was  ceasing  and  that  the  system  was  transi¬ 
tioning  from  outburst  to  quiescence.  This  is  also  illustrated 
by  Fig.  [lj  which  displays  the  X-ray  lightcurve  of  EXO  0748- 
676  as  observed  with  the  All-Sky  Monitor  (ASM)  onboard 
RXTE  since  1996.  The  decrease  in  source  activity  is  clearly 
seen  around  ~  4600  days. 

Chandra  observations  carried  out  in  2008  mid-October 
(i.e.,  after  the  transition  to  quiescence)  revealed  an  X-ray 
spectrum  composed  of  a  soft,  thermal  component  joined  by 
a  hard  powerlaw  tail  that  dominates  the  spectrum  above 
~  2  —  3  keV  (IDegenaar  et  al.  il2009l.  see  also  Sections  m  and 
E2]>.  This  is  frequently  seen  for  neutron  star  X-ray  bina- 


ries  in  quiescence  (e.f 

Rutledge  et  al.  19991;  in’t  Zand  et  al. 

200ll:  Tomsick  et  al. 

2004).  The  non-thermal  component  is 

usually  well-fitted  by  a  simple  powerlaw  with  index  1  2  (e.g., 
lAsai  et  all  119961') .  The  fractional  contribution  of  the  hard 
powerlaw  tail  to  the  0.5-10  keV  X-ray  flux  widely  varies 
amongst  sources  and  possibly  also  with  changing  luminosity 

d. Ionker  et,  alJEooil: 1.TonkeifeOOSlf .  The  physical  process  that 
yields  the  powerlaw  spectral  component  remains  elusive  (see 

e. g..  ICamnana  et  al.lll998l: ICampanall2003l'). 

Although  the  soft  spectral  component  has  been  ascribed 
to  low-level  accretion  onto  the  neutron  star  (IZampieri  et  ali 
1 199511 .  it  is  most  often  interpreted  as  thermal  surface  radia- 
tion  resulting  from  deep  crustal  heating  (iBrown  et  al.lfl998l). 
According  to  this  model,  the  accretion  of  matter  compresses 
the  neutron  star  crust,  which  induces  a  series  of  reactions 
that  deposit  heat  (e.g..lHaensel  fc  Zduniklli 99(ll.  [20031.  l200<-j: 
iGupta  et  al.ll2007h.  The  heat  energy  gained  via  this  process 
is  subsequently  lost  during  quiescent  episodes,  resulting  in  a 
candescent  thermal  emission  from  the  neutron  star  surface 
set  by  the  time-averaged  accretion  rate  of  the  system,  as 
well  as  the  rate  of  neutrino  emissions  from  the  stellar  core 
(e.g..  IBrown  et  aI1ll998l: |Colpi  et  al.ll200lh.  The  latter  cools 
the  interior  and  depends  on  the  equation  of  state  of  cold 
nuclear  matter  and  the  central  density  of  the  neutron  star 
(e.g..  lYakovlev  et  al.l l20().’j;  lYakovlev  fc  PethickH2004h.  Ob¬ 
servations  of  quiescent  neutron  stars  can  thus  provide  insight 
into  their  interior  properties. 

In  this  context,  a  special  opportunity  arises  for  the 
quasi-persistent  neutron  star  X-ray  binaries.  In  2001, 
both  KS  1731-260  and  MXB  1659-29  made  the  transi¬ 
tion  to  quiescence,  following  accretion  episodes  of  12.5 
and  2.5  years,  respectively  llWiinands  et  alJ  l200ll.  120021, 
120031.  120041;  ICackett  et  al.l  I2006I.  120081') .  More  recently,  in 
2007,  the  ~  1.6-year  long  outburst  of  XTE  J1701-462 
came  to  a  halt  (lAltamirano  et  al.ll2007l:  iHoman  et  al.  1I2007I: 

iFridriksson  et  al.ll2010l).  All  three  systems  were  subsequently 

monitored  with  Chandra  and  XMM-Newton,  which  revealed 
that  thermal  flux  and  neutron  star  temperature  were  gradu¬ 


Figure  1.  RXTE / ASM  20-day  averaged  lightcurve  (1.5-12  keV) 
of  EXO  0748-676,  illustrating  the  cessation  of  the  outburst  in  late 
2008.  The  dashed  vertical  line  corresponds  to  2008  December  31. 
The  arrows  indicate  the  times  of  our  four  sequences  of  Chandra 
observations,  which  were  performed  when  the  source  dropped  be¬ 
low  the  detection  limit  of  RXTE  (both  of  the  ASM  and  the  PCA). 

ally  decreasing  over  the  course  of  years  (see  also  Section  [4]). 
This  can  be  interpreted  as  cooling  of  the  neutron  star  crust 
that  has  been  heated  during  the  prolonged  accretion  out¬ 
burst  (IWiinands  et  aljfeooitlllutledge  et  al.ll2002l~).  Success¬ 
ful  modelling  of  the  observed  quiescent  X-ray  lightcurves 
with  neutron  star  thermal  evolution  models  supports  this 
hypothesis  and  provides  important  constraints  on  the  crust 
properties,  such  as  the  thermal  conductivity  llShternin  et  al.l 
120071;  iBrown  fe  Cummi~ngll2009l') . 

Along  these  lines  we  have  pursued  an  observational 
campaign  of  EXO  0748-676  to  study  the  time  evolution 
of  the  quiescent  X-ray  emission  following  its  long  accre¬ 
tion  outburst.  In  IDegenaar  et  al]  (l2009l~).  we  reported  on 
Chandra  and  Swift  observations  obtained  between  2008 
September  28  and  2009  January  30.  We  found  a  relatively 
hot  and  luminous  quiescent  system  with  a  temperature  of 
ZcT^j  ~  0.11  —  0.13  keV  and  a  thermal  0.01-100  keV  lumi¬ 
nosity  of  ~  (8  —  16)  x  1033  (d/7.4  kpc)2  erg  s_1.  No  clear 
decrease  in  effective  temperature  and  thermal  bolometric 
flux  was  found  over  the  five-month  time  span. 

I11  this  paper  we  report  on  continued  Swift  and  Chan¬ 
dra  observations  of  EXO  0748-676  during  its  quiescent  state. 
In  addition  we  include  an  archival  XMM-Newton  observa¬ 
tion  performed  ~  2  months  after  the  cessation  of  the  out- 
burst.  Previous  Chandra  and  Swift  observations  discussed 
bv  IDegenaar  et  all  (l2009h  were  re-analysed  in  this  paper  in 
order  to  obtain  a  homogenous  quiescent  lightcurve. 


2  OBSERVATIONS  AND  DATA  ANALYSIS 

Table  □  gives  an  overview  of  all  new  observations  of 
EXO  0748-676  discussed  in  this  paper.  A  list  of  earlier  Chan- 
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Table  1.  Observation  log. 


Satellite 

Obs  ID 

Date 

Exp.  time 
(ks) 

XMM 

0560180701* 

2008-11-06 

29.0  (MOS) 
22.9  (PN) 

Swift 

31272016 

2009-02-13 

3.5 

Swift 

31272017 

2009-02-20 

4.1 

Swift 

31272018* 

2009-02-23 

5.1 

Chandra 

9071* 

2009-02-23/24 

15.8 

10871* 

2009-02-25 

9.6 

Swift 

31272019* 

2009-03-01 

3.2 

Swift 

31272020 

2009-03-10 

5.1 

Swift 

31272021 

2009-03-16 

4.6 

Swift 

31272022 

2009-04-09 

3.5 

Swift 

31272023 

2009-04-16 

2.8 

Swift 

31272024 

2009-04-23 

4.8 

Swift 

31272025 

2009-05-07 

4.5 

Swift 

31272026 

2009-05-14 

3.6 

Swift 

31272027 

2009-05-28 

3.4 

Swift 

31272028 

2009-06-05 

4.1 

Chandra 

9072* 

2009-06-10 

27.2 

Swift 

31272029* 

2009-06-11 

4.3 

Swift 

31272030 

2009-06-18 

3.9 

Swift 

31272031* 

2009-06-26 

5.5 

Swift 

31272032* 

2009-07-03 

4.8 

Swift 

31272033 

2009-07-18 

5.5 

Swift 

31272034* 

2009-07-25 

5.8 

Swift 

31272035* 

2009-07-31 

10.3 

Swift 

31272036* 

2009-08-15 

9.4 

Swift 

31272037 

2009-08-25 

1.1 

Swift 

31272038 

2009-08-26 

7.4 

Swift 

31272039* 

2009-09-08 

4.7 

Swift 

31272040* 

2009-09-09 

4.3 

Swift 

31272041 

2009-10-01 

1.9 

Swift 

31272042 

2009-10-02 

1.8 

Swift 

31272043 

2009-10-07 

2.0 

Swift 

31272044 

2009-10-08 

2.4 

Swift 

31272045 

2009-10-09 

2.3 

Swift 

31272046* 

2009-11-05 

4.2 

Swift 

31272047* 

2009-12-21 

9.4 

Swift 

31272048* 

2010-10-01 

9.6 

Swift 

31272049 

2010-02-12/13 

11.3 

Swift 

31272050* 

2010-03-12/13 

9.5 

Chandra 

11059* 

2010-04-20 

27.4 

Note.  -  The 

observations 

marked  with  an 

asterisk  contain 

(part  of)  eclipses.  The  listed  exposure  times  represent  the 
duration  of  the  observations  uncorrected  for  eclipses. 

dr a  and  Swift-  observations  in  quiescence  can  be  found  in 
iDegenaar  et  all  (i2009i). 

2.1  XMM-Newton 

EXO  0748-676  was  observed  with  the  European  Photon 
Imaging  Camera  (EPIC)  onboard  XMM-Newton  on  2008 
November  6  from  08:30-16:42  ut  (see  also  iBassa  et  al.l 
120091).  Both  the  MOS  detectors  and  the  PN  camera  were 
operated  in  full  window  mode  and  using  the  medium  opti¬ 
cal  blocking  filter.  Data  reduction  and  analysis  was  carried 
out  with  the  Science  Analysis  Software  (SAS;  v.  9.0.0).  We 
reprocessed  the  Original  Data  Files  (ODF)  using  the  tasks 
EMPROC  and  EPPROC.  To  identify  possible  periods  of  high 
particle  background,  we  extracted  high-energy  lightcurves 


10  keV  for  the  MOS  and  between  10-12  keV  for  the 
PN).  No  strong  background  flares  occurred  during  the  ob¬ 
servation.  The  net  exposure  times  are  29.0  and  22.9  ks 
for  the  MOS  and  PN,  respectively.  EXO  0748-676  is  de¬ 
tected  at  count  rates  of  0.16  ±  0.01  counts  s_1  (MOS)  and 
0.55  ±  0.01  counts  s^1  (PN). 

Source  spectra  and  lightcurves  were  obtained  with  the 
software  task  evselect,  using  a  35  arcsec  circular  re¬ 
gion  and  applying  pattern  selections  0-12  and  0-4  for  the 
MOS  and  PN  data,  respectively.  Corresponding  background 
events  were  extracted  from  a  circular  region  with  a  radius 
70  arcsec.  For  the  MOS  cameras,  the  background  was  po¬ 
sitioned  on  a  source-free  region  on  the  same  CCD  as  the 
source.  For  the  PN  instrument  the  background  events  were 
extracted  from  an  adjacent  CCD,  at  the  same  distance  from 
the  readout  node  to  ensure  similar  low-energy  noise.  The  an¬ 
cillary  response  Hies  (arf)  and  redistribution  matrices  (rmf) 
were  generated  with  the  tasks  arfgen  and  rmfgen. 

The  EPIC  lightcurves  show  two  full  eclipses  (see  also 
IBassa  et  al  ]|20()iil ,  corresponding  to  eclipse  cycles  54384  and 
54385  in  the  numbering  system  oflParmar  et  al]  dl986l).  To 
calculate  the  correct  non-eclipse  time-averaged  fluxes,  we  re¬ 
duce  the  exposure  times  for  each  instrument  by  500  s  per 
eclipse,  which  is  the  approximate  length  of  the  eclipses  of 
EXO  0748-676  llParmar  et  al  .1  [Tfisfll:  I  Wolff  et  al.ll2009l).  Us¬ 
ing  the  tool  GRPPHA,  the  spectra  were  grouped  to  contain 
a  minimum  of  20  photons  per  bin.  The  combined  MOS  and 
PN  spectral  data  were  fitted  in  xspec  (v.  12)  in  the  0.5-10 
keV  energy  range,  with  all  spectral  parameters  tied  between 
the  three  EPIC  detectors. 

A  single  absorbed  powerlaw  model  provides  an  accept¬ 
able  fit  to  the  XMM-Newton  data  (red.  \2  =  1-3  for  466 
d.o.f.).  However,  the  spectral  index  is  unusually  large  for  an 
X-ray  binary  (P  =  4.7  ±  0.1)  and  suggests  that  the  spec¬ 
trum  has  a  thermal  shape.  A  simple  absorbed  blackbody 
model,  bbodyrad,  provides  an  adequate  fit  (x2  =  1-2  for 
466  d.o.f.),  although  the  inferred  emitting  region  has  a  much 
smaller  radius  than  expected  for  a  neutron  star  (~  2  —  4  km 
for  distances  of  5  —  10  kpc).  However,  it  is  thought  that  ra¬ 
diative  transfer  effects  in  the  neutron  star  atmosphere  cause 
the  emergent  spectrum  to  deviate  from  a  blackbody  (e.g., 
IZavlin  et  al.lll99^:  iR.utledge  et,  al.lll999l).  There  are  several 
neutron  star  atmosphere  models  available  within  XSPEC, 
which  yield  equivalent  results  (see  e.g..  iHeinke  et  aI1l200fil: 

I  Webb  fc  Barretll2007h.  In  the  remainder  of  this  work,  we  con¬ 
centrate  on  fitting  the  data  with  ameutron  star  atmosphere 
model  NSATMOS  (jHeinke  et  al.ll2006h. 

The  model  parameters  are  the  hydrogen  column  density 
(Ah),  the  neutron  star  mass  and  radius  (Mns  and  I?ns),  the 
effective  temperature  in  the  neutron  star  frame  (i.e.,  non- 
redshifted;  feT/ff),  the  source  distance  ( D )  and  a  normalisa¬ 
tion  factor,  which  parametrizes  the  fraction  of  the  surface 
that  is  radiating.  We  keep  the  latter  fixed  at  1  through¬ 
out  this  work,  which  corresponds  to  the  entire  neutron  star 
surface  emitting.  The  effective  temperature  as  seen  by  an 
observer  at  infinity  is  given  by  fcT/jf  =  kTe g  x  gr,  where 
gr  =  —  (2GMns)I(Nnsc2)  is  the  gravitational  redshift, 

G  the  gravitational  constant  and  c  the  speed  of  light. 

The  NSATMOS  model  provides  a  good  fit  to  the  data 
(X2  =  1.1  for  466  d.o.f.),  although  significant  residuals  above 
the  model  fit  are  present  for  energies  >2  —  3  keV.  We  model 
this  non-thermal  emission  by  adding  a  powerlaw  compo- 
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nent,  which  significantly  improves  the  fit  (\2  =  1.0  for  464 
d.o.f.;  an  F-test  suggests  a~  lx  10-14  probability  of  achiev¬ 
ing  this  level  of  improvement  by  chance).  Chandra  observa¬ 
tions  carried  out  in  2008  mid-October,  three  weeks  prior  to 
this  XMM-Newton  observation,  also  indicated  the  presence 
of  a  non-thermal  component  in  the  quiescent  spectrum  of 
EXO  0748-676  (IDegenaar  et  al.l  [20091').  Whereas  the  Chan¬ 
dra  data  could  not  constrain  the  powerlaw  index,  the  larger 
collective  area  of  XMM-Newton  provides  better  constraints 
for  the  fluxes  under  consideration. 

By  fitting  the  XMM-Newton  data  to  a  combined  NS  AT¬ 
MOS  and  powerlaw  model,  we  obtain  a  powerlaw  index  of 
r  =  1.7  ±  0.5,  i.e.,  in  between  the  values  of  T  =  1  and 
F  =  2  considered  by  IDegenaar  et  alJ  d2009l'l.  For  this  fit, 
which  furthermore  yields  1Vh  =  (7  ±  2)  x  1020  cnfi2  and 
_RNS  =  17.8  ±  1  km  for  fixed  values  of  D  =  7.4  kpc 
and  Mns  =  1.4  M@,  the  powerlaw  component  contributes 
~  10  percent  to  the  total  unabsorbed  0.5-10  keV  flux. 
This  is  lower  than  the  ~  15  —  20  percent  inferred  from 
the  Chandra  observations  performed  in  2008  mid-October 
(|Degenaar  et  al  am-  The  obtained  hydrogen  column  den¬ 
sity  is  consistent  with  values  found  for  EXO  0748-676  during 
its  outburst  (JVh  ~  (7  —  12)  x  1020  cm~2,  e.g..  ISidoli  et  al 
1200511 . 

2.2  Chandra 

We  obtained  three  new  Chandra/ ACIS-S  observations.  The 
first  consists  of  two  separate  exposures  obtained  on  2009 
February  23-24  22:07-03:15  UT  (obs  IP  9071)  and  2009 
February  25  12:32-15:59  UT  (obs  IP  10871),  lasting  for 
~  15.8  and  ~  9.6  ks,  respectively.  In  both  data  sets, 
EXO  0748-676  is  clearly  detected  at  a  count  rate  of  0.17  ± 
0.01  counts  s_1  (0.5-8  keV).  This  is  a  factor  ~  1.5  lower  than 
observed  in  2008  October,  when  the  source  was  detected 
with  Chandra /ACIS-S  at  a  rate  of  0.24  ±  0.01  counts  s~4. 
Two  full  eclipses  are  seen  in  the  liglitcurve  of  observation 
9071,  while  one  eclipse  is  present  in  that  of  10871  (eclipse 
cycle  numbers  55071,  55072  and  55080,  respectively). 

A  second  Chandra  observation  was  carried  out  on  2009 
June  10,  from  12:36-21:16  UT,  with  an  exposure  time  of 

27.2  ks  (obs  IP  9072).  In  this  observation  EXO  0748-676  is 
detected  at  a  count  rate  of  0.16  ±  0.01  counts  s_1  and  the 
lightcurve  shows  two  full  eclipses  (cycles  55740  and  55741). 
Furthermore,  a  27.4  ks  exposure  was  taken  on  2010  April  20 
from  02:37-11:28  UT  (obs  IP  11059),  which  captured  three 
full  eclipses  (see  Fig.  [2]  these  correspond  to  eclipse  cycle 
numbers  57708,  57709  and  57710),  and  detected  the  source 
at  a  count  rate  of  0.14  ±  0.01  counts  s_1.  Similar  to  our 
treatment  of  the  XMM-Newton  data,  we  reduce  the  expo¬ 
sure  times  of  all  Chandra  observations  by  500  s  per  eclipse. 
There  are  no  indications  of  background  flares,  so  the  full 
data  set  is  used  in  further  analysis. 

We  reduced  the  data  employing  the  CIAO  tools  (v.  4.2) 
and  standard  Chandra  analysis  threads.  For  all  three  ob¬ 
servation  sequences,  the  ACIS-S3  CCP  was  operated  in  a 
1/8  sub-array,  resulting  in  a  frame-time  of  0.4  s.  For  the  ob¬ 
served  count  rates,  the  pile-up  fraction  is  <  2  percent,  so 
no  further  corrections  were  made.  Spectra  were  extracted 
with  the  tool  psextract,  while  the  rmf  and  arf  files  were 
created  using  mkacisrmf  and  mkarf,  respectively.  We  em¬ 
ploy  a  circular  region  of  3  arcsec  to  obtain  source  events  and 


Figure  2.  Chandra/ ACIS-S  0.5—8  keV  lightcurve  of  EXO  0748- 
676  obtained  on  2010  April  20  (obs  ID  11059).  Each  point  repre¬ 
sents  200  s  of  data.  Three  eclipses  are  visible. 

a  10-25  arcsec  annulus  for  the  background.  We  also  repro- 
cessed  the  Chandra  observation  obtained  in  2008  October 
(see  IPegenaar  et  alJ  120091')  to  benefit  from  the  calibration 
update  that  was  released  in  2009  Pecember.  Prior  to  spec¬ 
tral  fitting,  the  spectra  were  grouped  to  contain  at  least  20 
photons  per  bin.  We  find  no  spectral  differences  between 
the  two  separate  exposures  performed  in  2009  February  and 
therefore  we  tie  all  parameters  for  these  two  spectra  in  the 
main  fits. 

The  Chandra  observations  obtained  in  2009  February 
and  June  are  well-fitted  by  a  single-component  thermal 
model,  but  the  2010  April  data  shows  significant  residuals 
at  energies  >2  —  3  keV.  If  we  include  a  powerlaw  with  pho¬ 
ton  index  F  =  1.7,  as  found  from  fitting  the  XMM-Newton 
data  (see  Section  12.111.  this  model  component  contributes 
~  10,  ~  5  and  ~  15  percent  to  the  total  unabsorbed  0.5- 
10  keV  flux  for  the  data  taken  in  2009  February,  June  and 
2010  April,  respectively.  Fig.  [3] compares  the  Chandra  spec¬ 
tral  data  obtained  on  2008  October  and  2010  April,  showing 
that  both  spectral  components  decreased  over  the  18-month 
time  span  that  separates  the  two  observations. 

2.3  Swift 

In  addition  to  the  Chandra  and  XMM-Newton  observations, 
we  have  been  monitoring  EXO  0748-676  on  a  regular  basis 
with  the  XRT  aboard  the  Swift  satellite.  Starting  in  2008 
late-September,  approximately  2  —  3  pointings  were  per¬ 
formed  each  month  with  a  typical  duration  of  ~  3  —  5  ks 
per  observation,  and  a  separation  of  ~  1  —  2  weeks.  From 
2009  November  onwards,  the  cadence  was  lowered  to  one  ob¬ 
servation  per  month  with  a  longer  exposure  time  when  pos¬ 
sible  (see  Table  [l]).  EXO  0748-676  is  detected  in  the  XRT 
observations  at  count  rates  of  ~  (1  —  5)  x  10~2  counts  s_1. 

All  Swift /XRT  observations  were  obtained  in  the 
photon-counting  (pc)  mode  and  were  processed  using  the 
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Figure  3.  Spectra  of  the  Chandra  observations  of  2008  October 
(black)  and  2010  April  (red),  along  with  the  model  fits  (solid 
lines).  The  dotted  and  dashed  lines  represent  the  contributions 
of  the  NSATMOS  and  powerlaw  components  for  the  2008  and  2010 
data,  respectively. 

xrtpipeline  with  standard  quality  cuts  (event  grade  0-12). 
Using  xselect  (v.  2.4),  we  extracted  source  spectra  from 
a  circular  region  with  a  radius  of  35  arcsec  (~  15  pixels), 
which  optimises  the  signal  to  noise  ratio  at  the  observed 
count  rates  fevans  et  all  12007).  Corresponding  background 
events  were  averaged  over  three  source-free  regions  of  sim¬ 
ilar  shape  and  size.  Employing  the  tool  xrtexpomap,  we 
created  exposure  maps  to  account  for  the  effective  area  of 
the  CDD,  while  arfs  generated  with  xrtmkarf  account  for 
vignetting  and  point-spread-function  corrections.  The  latest 
rmf  (v.  11)  was  obtained  from  the  CALDB  database. 

We  reduced  the  exposure  times  of  those  observa¬ 
tions  that  contained  eclipses  according  to  the  ephemeris  of 
IWolff  et  ah  (2009)  to  calculate  the  correct  non-eclipse  time- 
averaged  fluxes  (see  Table  [lj .  Furthermore,  Swift  observa¬ 
tions  obtained  within  a  2-day  time  span  were  summed  to 
improve  the  data  statistic^.  This  seems  justified,  since  the 
Chandra  data  do  not  reveal  any  spectral  changes  on  such 
time  scales  (see  Section  T2. 2 1) . 

3  RESULTS 
3.1  Spectral  fits 

We  first  fit  the  Chandra  and  XMM-Newton  data  simulta¬ 
neously  to  an  absorbed  NSATMOS  model  joined  by  a  power- 
law,  to  explore  the  best-fit  values  for  the  neutron  star  mass 
and  radius,  source  distance  and  hydrogen  column  density. 
We  include  the  first  set  of  Chandra  observations  obtained 
in  2008  October  (discussed  in  iDegenaar  et,  ahl  120091)  in  the 
analysis.  For  the  neutral  hydrogen  absorption,  we  use  the 
PHABS  model  with  the  default  xspec  abundances  and  cross- 
sections.  The  powerlaw  index  is  fixed  to  T  =  1.7  (the  best 
fit- value  obtained  from  XMM-Newton  observations;  see  Sec¬ 
tion  12.11) ,  because  there  are  not  sufficient  counts  at  higher 

1  This  is  the  case  for  obs  IDs  31272037/38,  31272039/40  and 
31272043/44/45;  see  Table [T] 


energies  in  either  the  Swift  or  the  Chandra  spectra  to  allow 
this  component  to  vary.  The  powerlaw  normalisation  is  left 
as  a  free  parameter. 

If  the  neutron  star  mass  and  radius  are  fixed  to  canon¬ 
ical  values  of  Mns  =  1.4  Mq  and  R  =  10  km,  and 
in  addition  the  source  distance  is  fixed  to  the  best  esti- 
mate  obtained  from  type-I  X-ray  bursts  ( D  =  7.4  kpc; 
iGallowav  et  al.1 120081) ,  the  hydrogen  column  density  pegs 
at  its  lower  limit  (Ah  =  0).  When  the  distance  is  left  to 
vary  freely,  the  best-fit  value  is  4.6  ±  0.3  kpc,  which  is  just 
outside  the  range  obtained  from  X-ray  burst  analysis  (5- 
8.3  kpc:IOallowav  et  al.ll2008).  Therefore,  we  choose  to  keep 
the  distance  fixed  at  7.4  kpc,  and  instead  allow  the  neutron 
star  radius  to  vary.  This  way,  we  obtain  best-fit  values  of 
Ah  =  (7±1)  x  1020  cm-2  and  R  =  15.6±0.8  km.  If  addition¬ 
ally  the  neutron  star  mass  is  left  free  to  vary  in  the  fit,  this 
parameter  is  not  strongly  constrained  (Mns  ~  1.6±0.6  Mq). 
In  the  final  fits  we  choose  to  fix  the  neutron  star  mass  to 
Mns  =  1.4  Mq,  because  otherwise  the  uncertainty  in  this 
quantity  will  dominate  the  errors  of  the  other  parameters. 

The  Swift  data  do  not  provide  sufficient  statistics  to 
constrain  the  powerlaw  component.  We  do  include  a  power- 
law  in  the  fits,  but  fix  both  the  index  and  the  normalisation 
of  this  component.  Since  it  is  unclear  how  the  powerlaw 
exactly  evolves  over  time,  we  adjust  the  powerlaw  normal¬ 
isation  for  the  Swift  observations  such  that  it  always  con¬ 
tributes  10  percent  of  the  total  unabsorbed  0.5-10  keV  flux. 
After  treating  each  Swift  observation  separately,  we  found 
that  the  thermal  flux  and  neutron  star  temperature  did  not 
evolve  significantly  between  consecutive  observations.  To  im¬ 
prove  the  statistics,  we  therefore  sum  the  Swift  data  into 
groups  spanning  ~  1  —  4  weeks  of  observations,  resulting  in 
exposure  times  of  ~  10  —  20  ks  (see  Table  0.  These  spectra 
were  grouped  to  contain  a  minimum  of  20  photons  per  bin. 

For  the  final  spectral  analysis,  we  fit  all  XMM-Newton, 
Chandra  and  Swift  data  with  an  absorbed  NSATMOS  plus 
powerlaw  model,  where  Ah  =  7xl020  cm-2,  Mns  =  1.4  Mq, 
J?ns  =  15.6  km,  D  =  7.4  kpc  and  T  =  1.7  are  fixed,  while 
the  neutron  star  effective  temperature  is  left  as  a  free  pa¬ 
rameter.  The  powerlaw  normalisation  is  fixed  for  the  Swift 
data,  but  left  to  vary  freely  for  the  Chandra  and  XMM- 
Newton  observations.  We  fit  all  data  in  the  0.5-10  keV  en¬ 
ergy  range  and  deduce  the  absorbed  and  unabsorbed  fluxes 
in  this  band.  The  thermal  model  fit  is  extrapolated  to  the 
energy  range  of  0.01-100  keV  to  estimate  the  thermal  bolo- 
metric  flux.  The  results  from  fitting  the  X-ray  spectra  in  this 
way  are  presented  in  Table0  The  effective  temperatures  and 
thermal  bolometric  fluxes  derived  from  Chandra,  Swift  and 
XMM-Newton  data  are  displayed  in  Fig.  [4]  Examination  of 
Fig.  0  suggests  that  there  is  a  small  but  discernible  offset  in 
the  thermal  flux  and  neutron  star  temperature  as  deduced 
from  the  different  satellites.  This  is  briefly  discussed  in  Sec¬ 
tion  [373] 

3.2  Lightcurve  fits 

Fig.  [4]  clearly  reveals  a  decaying  trend  in  thermal  flux  and 
temperature.  To  investigate  the  decay  shape,  we  fit  the  tem¬ 
perature  curve  with  an  exponential  decay  function  of  the 
form  y(t)  =  a  e_v  *°)/r,  where  a  is  a  normalisation  con¬ 
stant,  to  is  the  start  time  of  the  cooling  curve  and  r  the 
e-folding  time.  Given  the  offset  between  the  different  in- 
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Table  2.  Results  from  fitting  the  spectral  data. 


Satellite 

Date 

At 

(days) 

Pow.  frac. 
(%) 

kT% 

(eV) 

Fx 

/?th 

^bol 

-^bol 

Red.  x2 
(d.o.f.) 

Swift f 

2008-09-28  - 

2008-10-07 

4.9 

10  fix 

123.7  ±5.4 

1.31  ±0.22 

1.53  ±0.26 

10.0  ±  1.7 

0.93  (8) 

Chandra f 

2008-10-12/13/15 

1.4 

20  ±3 

118.8  ±  0.9 

1.23  ±0.02 

1.31  ±  0.04 

8.6  ±0.3 

1.03  (175) 

Swift] 

2008-10-29  - 

2008-11-02 

2.2 

10  fix 

118.3  ±2.6 

1.10  ±0.09 

1.28  ±0.11 

8.4  ±0.7 

0.67  (14) 

XMM 

2008-11-06 

0.2 

7  ±  2 

120.7  ±  0.4 

1.14  ±0.01 

1.39  ±0.02 

9.1  ±0.1 

1.08  (467) 

Swift] 

2008-11-28  - 

2008-12-20 

11.0 

10  fix 

118.7  ±2.6 

1.11  ±0.10 

1.30  ±0.12 

8.5  ±0.8 

1.09  (14) 

Swift f 

2009-01-10  - 

2009-01-30 

9.8 

10  fix 

116.2  ±2.2 

0.99  ±0.07 

1.19  ±0.09 

7.8  ±0.6 

1.09  (18) 

Swift 

2009-02-13  - 

2009-02-23 

5.1 

10  fix 

117.2  ±2.4 

1.02  ±0.08 

1.23  ±0.10 

8.1  ±0.7 

0.50  (15) 

Chandra 

2009-02-23/25 

0.9 

12  ±4 

113.5  ±  1.3 

0.91  ±0.03 

1.09  ±0.05 

7.1  ±0.3 

0.90  (139) 

Swift 

2009-03-01  - 

2009-03-16 

7.2 

10  fix 

115.6  ±2.3 

0.97  ±0.08 

1.17  ±  0.09 

7.7  ±0.6 

0.79  (16) 

Swift 

2009-04-09  - 

2009-04-23 

7.1 

10  fix 

112.2  ±  2.8 

0.86  ±0.09 

1.03  ±0.11 

6.8  ±0.7 

1.16  (10) 

Swift 

2009-05-07  - 

2009-06-05 

14.7 

10  fix 

114.2  ±2.3 

0.92  ±0.07 

1.11  ±  0.09 

7.3  ±0.6 

0.88  (16) 

Chandra 

2009-06-10 

0.2 

4  ±  3 

111.0  ±0.7 

0.75  ±0.01 

0.99  ±  0.03 

6.5  ±0.2 

1.19  (93) 

Swift 

2009-06-11  - 

2009-07-03 

11.5 

10  fix 

111.9  ±2.2 

0.75  ±0.07 

1.03  ±0.08 

6.7  ±0.5 

0.99  (17) 

Swift 

2009-07-18  - 

2009-07-31 

6.9 

10  fix 

110.5  ±  2.1 

0.79  ±0.06 

0.98  ±0.07 

6.4  ±0.5 

0.42  (18) 

Swift 

2009-08-15  - 

2009-09-09 

12.7 

10  fix 

110.0  ±  1.8 

0.78  ±0.05 

0.96  ±  0.06 

6.3  ±0.4 

1.30  (26) 

Swift 

2009-10-01  - 

2009-11-05 

17.4 

10  fix 

108.0  ±2.6 

0.69  ±0.07 

0.88  ±  0.09 

5.8  ±0.6 

1.49  (11) 

Swift 

2009-12-21  - 

2010-10-01 

10.2 

10  fix 

109.4  ±2.0 

0.74  ±0.06 

0.94  ±0.07 

6.1  ±0.5 

1.51  (19) 

Swift 

2010-02-12  - 

2010-03-13 

15.0 

10  fix 

109.4  ±2.0 

0.76  ±0.06 

0.94  ±0.07 

6.1  ±0.5 

1.25  (20) 

Chandra 

2010-04-20 

0.2 

15  ±4 

108.6  ±  1.1 

0.77  ±0.02 

0.91  ±  0.04 

6.0  ±0.2 

0.79  (91) 

Note.  -  The  observations  marked  by  a  dagger  were  already  discussed  in  IDegenaar  et  all  f 200ft) .  but  re-fitted  in  this  work.  These 
results  were  obtained  by  using  a  combined  absorbed  NSATMOS  and  powerlaw  model,  where  ,Yh  =  7  X  1020  cm-2,  Mxs  =  1.4  Mq, 
/i\S  =  15-6  km,  D  =  7.4  kpc  and  T  =  1.7  were  kept  fixed.  The  quoted  errors  represent  90%  confidence  levels.  Fx  represents  the 
0.5-10  keV  total  model  flux  and  Ff^ j  gives  the  0.01-100  keV  NSATMOS  flux;  both  are  unabsorbed  and  in  units  of  10-12  erg  cm-2  s-1. 
Lbol  gives  the  0.01-100  keV  luminosity  of  the  NSATMOS  model  component  in  units  of  1033  erg  s~  ]  and  assuming  a  source  distance 
of  D=7.4  kpc.  At  represents  the  time  interval  of  the  observations  in  days  and  the  fractional  powerlaw  contribution  is  given  in  a 
percentage  of  the  total  unabsorbed  0.5-10  keV  flux. 


struments  (see  Section  [3.31).  we  perform  different  fits  to  the 
Chandra  and  Swift  data.  We  fix  to  to  2009  September  5 
(MJD  54714),  which  is  in  between  the  first  non-detection  by 
RXTE /PCA  and  the  first  Swift  observation  of  the  source 
(see  IDegenaar  et  alJ 120091) . 

This  simple  exponential  decay,  represented  by  the  dot¬ 
ted  lines  in  Fig.  [5]  yields  an  e-folding  time  of  6121.7  ± 
2004.0  days  for  the  Chandra  data,  but  does  not  provide  a 
good  fit  (red.  \2  —  6.0  for  2  d.o.f.).  For  the  Swift  lightcurve 
we  find  r  =  5328.1  ±674.7  days  (red.  ,\2  =  0.5  for  12  d.o.f.). 
If  we  include  a  constant  offset  (i.e.,  y(t)  =  a  e-0-*o)/'r  -|-5)) 
we  obtain  a  better  fit  for  the  Chandra  data  (solid  lines  in 
Fig.  0,  yielding  a  normalisation  of  13.4  ±  0.2  eV,  an  e- 
folding  decay  time  of  191.6  ±  9.7  days  and  a  constant  level 
of  107.9  ±  0.2  eV  (red.  y2  =  0.02  for  1  d.o.f.).  For  the  Swift 
data  this  results  in  a  =  17.2  ±  1.8  eV,  r  =  265.6  ±  100.0  days 
and  b  =  106.2  ±  2.5  eV  (red.  \2  —  0.34  for  11  d.o.f.),  which 
is  consistent  with  the  Chandra  fit. 

Although  an  exponential  decay  provides  an  adequate 
description  of  the  data  of  EXO  0748-676,  as  has  been  found 
for  other  sources  (e.g.JCackett  et,  al.  2006;  Fridriksson  et  all 
l2010h.  mathematically  a  neutron  star  crust  is  expected 
to  cool  via  a  (broken)  powerlaw  llEichler  fe  Chenall989l: 
iBrown  fc  Cummingll2009ri.  If  we  fit  a  single  powerlaw  of  the 
form  y(t )  =  A(t—to)B  to  the  Chandra  data,  we  find  an  index 
of  —0.03  ±  0.01  and  a  normalisation  of  134.4  ±  1.0  eV  (red. 

=  0.13  for  2  d.o.f.).  For  the  Swift  observations  we  find 
B  =  -0.05  ±0.01  and  A  =  144.7  ±3.8  eV  (red.  *2  =  0.4  for 
12  d.o.f.).  These  powerlaw  fits  are  indicated  by  the  dashed 
lines  in  Fig.  [5] 

A  broken  powerlaw  also  yields  an  acceptable  fit  to  the 


Swift  data  (red.  %2  =  0.3  for  10  d.o.f.).  We  find  a  normali¬ 
sation  of  A  =  135.0  ±  17.8  eV,  a  break  at  166.0  ±  99.2  days 
and  decay  indices  of  —0.03  ±  0.03  and  —0.06  ±  0.02  before 
and  after  the  break,  respectively.  This  fit  is  indicated  by 
the  dashed-dotted  curve  in  Fig.  [5]  There  are  not  sufficient 
Chandra  observations  to  fit  a  broken  powerlaw  decay.  We 
note  that  the  shape  of  the  decay  curve  of  EXO  0748-676 
is  not  strongly  affected  by  our  choice  of  spectral  parame¬ 
ters  (Ah,  Mns,  72ns,  and  T)  or  assumed  distance  (see  also 
previous  studies  by  e.g..lWiinands  et  al.ll2004l;ICackett  et  al.l 

l2008h. 


3.3  Instrument  cross-calibration 

The  quiescent  lightcurve  presented  in  Fig.  [4]  indicates  that 
the  thermal  flux  and  temperature  inferred  from  the  XMM- 
Newton  observation  lie  exactly  on  the  trend  indicated  by 
the  Swift  data  points.  However,  the  Chandra  results  are  sys¬ 
tematically  lower.  This  apparent  offset  (~  6  percent  for  the 
flux  lightcurve)  is  possibly  due  to  cross-calibration  issues  be¬ 
tween  the  two  satellites  (see  e.g..lKirsch  et  al1l2005l).  Such  a 
discrepancy  has  not  been  found  for  Swift  and  XMM-Newton 
(e.g.,  IKirsch  et  al.l  [2005") .  which  appears  to  be  reflected  in 
our  results  as  well.  We  note  that  regardless  of  the  offset  be¬ 
tween  the  Chandra  and  Swift  points,  both  data  sets  show  a 
similar  decay  in  thermal  flux  and  neutron  star  temperature 
(see  Section  l±2li . 
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Figure  5.  Evolution  of  the  effective  temperature  of  EXO  0748-676  fitted  to  different  decay  functions  (see  Section  m-  The  left  image 
displays  Chandra  data  and  exponential  decay  fits  both  with  and  without  a  constant  offset  (solid  and  dotted  line,  respectively),  as  well 
as  a  decaying  powerlaw  (dashed  curve).  The  right  image  shows  Swift  observations,  where  the  dashed  line  is  again  a  powerlaw  fit,  while 
the  solid  and  dotted  curves  are  exponential  decays.  In  addition,  this  plot  includes  a  fit  to  a  broken  powerlaw,  which  is  represented  by 
the  dashed-dotted  line. 
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Figure  4.  Evolution  of  the  bolometric  flux  (top)  and  effective 
temperature  (bottom)  of  EXO  0748—676,  deduced  from  Chan¬ 
dra  /ACIS-S  (black  squares),  Swift /X RT  (grey  triangles)  and 
X M M-  Ne. w t a n /EPIC  (black  star)  data.  Multiple  Swift  obser¬ 
vations  were  summed  to  improve  the  data  statistics  (see  Sec¬ 
tion  12.311 . 


4  DISCUSSION 

We  discuss  Chandra,  Swift  and  XMM-Newton  observations 
obtained  after  the  cessation  of  the  very  long  (~24  year) 
active  period  of  EXO  0748-676.  The  data  cover  the  first 
19  months  (1.6  years)  of  the  quiescent  phase.  By  assum¬ 


ing  a  neutron  star  atmosphere  model  NSATMOS,  we  find  a 
clear  decrease  in  neutron  star  effective  temperature  from 
k.Tftf  ~  124  to  109  eV.  The  thermal  bolometric  flux  decayed 
from  Fbo!  ~  1.5  x  10-12  to  0.9  x  10-12  erg  cm-2  s-1. 


In  addition  to  a  soft,  thermal  component,  the  Chandra 
and  XMM-Newton  observations  show  evidence  for  a  hard 
powerlaw  tail  with  index  T  =  1.7.  The  fractional  contribu¬ 
tion  of  the  hard  spectral  component  to  the  total  unabsorbed 
0.5-10  keV  flux  initially  decreased  from  ~  20  percent  in 
2008  October  to  ~  4  percent  in  2009  June.  However  obser¬ 
vations  carried  out  in  2010  April  suggest  that  the  powerlaw 
fraction  increased  again  to  ~  15  percent.  Similar  behav¬ 
ior  has  been  observed  for  several  other  quiescent  neutron 
star  systems  lljonker  et  al.ll2003:IJonker|2008t),  although  oth- 
ers  show  more  irregular  behaviour  (Fridriksson  et  al.  120101). 
In  Cen  X-4,  the  powerlaw  tail  in  the  quiescent  spectrum 
shows  variations  that  appear  to  be  linked  to  changes  in  the 
thermal  component,  possibly  caused  by  low-level  accretion 
llCackett  et  al.ll201Cl). 


The  gradual  decrease  in  thermal  flux  and  neutron  star 
temperature  observed  for  EXO  0748-676  can  be  interpreted 
as  the  neutron  star  crust  cooling  down  in  quiescence  af¬ 
ter  it  has  been  heated  during  its  long  accretion  outburst. 
Fig.  [6]  compares  our  data  of  EXO  0748-676  with  the  crust 
cooling  curves  observed  for  the  neutron  star  X-ray  binaries 
KS  1731-260,  MXB  1659-29  and  XTE  J1701-462.  This  plot 
clearly  shows  that  the  quiescent  lightcurve  of  EXO  0748-676 
is  relatively  flat  compared  to  the  others.  We  have  observed 
our  target  over  the  first  ~  19  months  after  the  cessation 
of  the  outburst  and  during  this  time  the  thermal  bolomet¬ 
ric  flux  has  decreased  by  a  factor  of  ~  1.7.  In  a  similar 
time  span,  the  thermal  bolometric  fluxes  of  KS  1731-260, 
MXB  1659-29  and  XTE  J1701-462  had  decreased  by  a  fac¬ 
tor  of  ~  3.5,  6  and  2.5,  respectively  (seelCackett  et  ahll2006l: 
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lEridriksson  et  all  |2()I  (Jl .  The  effective  neutron  star  temper¬ 
ature  of  EXO  0748-676  has  decreased  by  about  10%,  com¬ 
pared  to  ~  30,  40  and  20%  for  KS  1731-260,  MXB  1659-29 
and  XTE  J1701-462. 

Although  the  observed  fractional  changes  in  neutron 
star  temperature  and  thermal  bolometric  flux  are  smaller 
for  EXO  0748-676  than  for  the  other  three  sources,  the 
decay  itself  may  not  be  markedly  different.  The  quiescent 
lightcurves  of  KS  1731-260,  MXB  1659-29  and  XTE  J1701- 
462  can  be  fit  with  an  exponential  decay  function  level¬ 
ling  off  to  a  constant  value,  yielding  e-folding  times  of 
~  305  ±  50,  ~  465  ±  25  and  ~  120  ±  25  days,  respectively 
llCackett  et  al.ll2008l:  iFridriksson  et  al.l[2010h .  For  the  Chan¬ 
dra  data  of  EXO  0748-676,  we  find  an  e-folding  time  of 
~  192  ±10  days  (see  Section  [372]) .  These  decay  times  provide 
a  measure  of  the  thermal  relaxation  time  of  the  neutron  star 
crust,  which  depends  on  the  composition  and  structure  of 


ness  of  the  crust  (e. 

i..  Lattimer  et  al.  1994!;  Rutledge  et  al. 

2002|;  Shternin  et  al. 

20Q21:  Brown  &  Cumming  2009). 

culate  theoretical  cooling  curves  for  KS  1731-  260,  assum¬ 
ing  different  physics  for  the  crust  and  core.  These  authors 
present  simulations  for  both  an  amorphous  crust  and  an 
ordered  crystalline  lattice.  For  the  latter,  the  spread  of 
nuclide  charge  numbers  (Z)  in  the  crust  matter  is  small, 
which  is  referred  to  as  a  low  level  of  impurities,  and  re¬ 
sults  in  a  highly  conductive  crust.  A  large  number  of  im¬ 
purities  gives  an  amorphous  structure,  which  affects  the 
thermal  properties  of  the  crust  and  results  in  a  low  con¬ 
ductivity.  In  addition,  IR.utledge  et  al„  (20021  explore  stan¬ 
dard  (i.e.,  slow)  and  enhanced  neutrino  cooling  mech¬ 
anisms,  yielding  different  core  temperatures.  Comparing 
our  results  on  EXO  0748-676  with  the  decay  shapes  re¬ 
sulting  from  those  calculations  suggests  that  the  neutron 
star  has  a  highly  conductive  crust,  similar  to  what  has 
been  inferred  for  the  other  three  sources  |Wiinandsjrt_alJ 


2001  12004  ICaekett  et  alJ  120061:  IShternin  et  al.1  120071: 
Brown  fe  Cummindl2009r  Fridriksson  et  al.l  2Qlol).  The  fact 
that  the  decay  curve  of  EXO  0748-676  is  rather  flat  may  be 
explained  in  terms  of  a  relatively  small  temperature  gradi¬ 
ent  and  thus  lower  thermal  flux  across  the  core-crust  bound- 
ary  (compare  the  model  curves  of  IR.utledge  et  ai1l2002l.  for 
a  highly  conductive  crust  and  different  core  temperatures). 
This  can  be  due  to  a  combination  of  a  warm  neutron  star 
core  (see  below)  and  a  relatively  low  mass-accretion  rate 
during  outburst. 

The  exponential  decay  fit  to  the  Chandra  data  of 
EXO  0748-676  indicates  that  the  neutron  star  crust  might 
already  be  close  to  restoring  equilibrium  with  the  core.  The 
fit  results  in  a  quiescent  base  level  of  107.9  ±  0.2  eV,  while 
we  found  a  temperature  of  108.6  ±  1.1  eV  for  the  obser¬ 
vation  performed  in  2010  April.  Prior  to  its  last  outburst, 
EXO  0748-676  was  observed  in  quiescence  with  the  EIN¬ 
STEIN  observatory,  displaying  a  0.5-10  keV  unabsorbed 
flux  of  (8. 4+17)  x  10-13  erg  cm-2  s_1  jGarcia  fe  Callananl 
1999).  Our  Chandra  observations  of  2010  April  detected 
EXO  0748-676  at  a  0.5-10  keV  unabsorbed  flux  of  (7.7  ± 
0.2)  x  10~13  erg  cm-2  s_1  (see  Tablc[2]|.  Assuming  that  the 
EINSTEIN  detection  caught  EXO  0748-676  at  its  quies¬ 
cent  base  level,  this  supports  the  idea  that  the  crust  has 
nearly  cooled  down.  This  would  imply  that  the  neutron  star 


in  EXO  0748-676  is  relatively  hot  Ccf.  Iffeinke  et  al]l2009h. 
suggesting  that  either  standard  cooling  mechanisms  are  op¬ 
erating  in  the  core  and  that  the  neutron  star  is  not  very 
massive,  or  that  the  time-averaged  mass- accretion  rate  of 
the  system  is  very  high. 

The  nuclear  reactions  induced  in  the  crust  of  an  ac¬ 
creting  neutron  star  are  expected  to  generate  a  thermal  lu- 
minosity  of  Lbol  ~  ( M)Qnuc/mu ,  once  the  accretion  halts 
fe.g..  [Brown  et  all  Il998l:  IColpi  et  all  1200111.  Here,  Qnuc  ~ 
2  MeV  is  the  nuclear  energy  deposited  per  accreted  baryon 
dGupta  et  al.l  120071:  iHaensel  &i  Zdunikl  l2008h.  mu  is  the 
atomic  mass  unit  and  (M)  is  the  time-averaged  accretion 
rate  of  the  system.  The  latter  can  be  expressed  as  (M)  = 
(. Mob )  x  tob/irec,  where  {Mob)  is  the  average  accretion  rate 
during  outburst  episodes,  t0b  is  the  outburst  duration  and 
tree  is  the  system’s  recurrence  time.  The  factor  t0\,/tTec  rep¬ 
resents  the  duty  cycle  of  the  system. 

The  outburst  of  EXO  0748-676  started  between  1980 
May  and  1984  July  (see  Section  [T])  and  the  system  re¬ 
turned  to  quiescence  in  2008  September,  i.e.,  t0b  =  24  —  28 
years.  The  estimated  average  bolometric  flux  during  out- 
burst  is  ~  (7  —  8)  x  1Q~10  erg  cm-2  s_1  dSidoli  et  al.l 
120051:  iBoirin  et  al.l  120071).  For  an  accretion  luminosity  of 
Ls.cc  =  (G,MNS/-RNs)(M0b},  this  translates  into  a  mass- 
accretion  rate  of  (Mob)  ~  3  x  1016  g  s^1  for  a  canoni¬ 
cal  neutron  star  with  M  =  1.4  Mq  and  R  =  10  km  (as¬ 
suming  a  distance  of  D  =  7.4  kpc)0  As  argued  above, 
EXO  0748-676  might  already  be  close  to  reaching  its  quies¬ 
cent  base  level.  To  explain  a  quiescent  bolometric  luminos¬ 
ity  of  ~  6  x  1033  (D/7.4  kpc)2  erg  s_1  (as  measured  during 
the  Chandra  observation  of  2010  April)  in  terms  of  thermal 
emission  resulting  from  deep  crustal  heating  (i.e.,  opposed 
to  continued  accretion),  it  requires  that  EXO  0748-676  has 
a  duty  cycle  of  ~  10  percent.  This  would  imply  an  expected 
recurrence  time  of  ~  200  —  250  years  if  the  observed  outburst 
is  typical  for  the  long-term  behaviour  of  this  source. 


_ Brownet^L  dl99Sl~).  iRutledge  et  al.l  d2000l'l  and 

IColpi  et  al.l  I200H)  have  suggested  that  EXO  0748-676 


continues  to  accrete  in  quiescence,  because  the  quiescent 
luminosity  inferred  from  the  1980  EINSTEIN  observation  is 
higher  than  predicted  by  standard  cooling  models.  However, 
these  conclusions  are  based  on  an  assumed  duty  cycle  of 
~  1  percent,  but  we  have  no  a  priori  knowledge  about 
this.  Although  we  cannot  exclude  that  the  system  is  indeed 
accreting  in  quiescence,  the  above  estimates  show  that  a 
duty  cycle  of  ~  10  percent  (which  is  not  uncommon  for 
neutron  star  transients)  can  explain  the  observed  quiescent 
level  of  EXO  0748-676  as  being  due  to  thermal  emission 

resulting  from  deep  crustal  heating.  _ 

Recently,  iBrown  fe  Cummind  l|2009ll  demonstrated 
that  the  cooling  of  a  neutron  star  crust  is  expected  to  fol¬ 
low  a  broken  powerlaw  decay.  A  break  is  predicted  to  occur 
due  to  a  transition  in  the  crystal  structure  of  the  crust  mat¬ 
ter,  and  the  slope  before  the  break  reflects  the  heat  flux  from 
the  outer  crustal  layers.  Therefore,  we  also  fitted  the  neutron 


2  We  note  that  EXO  0748—676  is  an  eclipsing  system  and  there¬ 
fore  part  of  the  central  X-ray  flux  may  be  intercepted  from  our 
line  of  sight.  However,  the  X-ray  burst  behaviour  of  the  source  is 
consistent  with  the  mass-accretion  rate  inferred  from  the  observed 
X-ray  luminosity  llBoirin  et  al.  I|2007|). 
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star  temperatures  obtained  for  EXO  0748-676  to  a  power- 
law  and  found  decay  indices  of  — 0.03±0.01  and  — 0.05±0.01 
for  the  Chandra  and  Swift  data  sets,  respectively.  The  Swift 
observations  indicate  that  a  possible  break  in  the  quiescent 
lightcurve  may  have  occurred  ~  67  —  265  days  after  the  ces¬ 
sation  of  the  outburst  (see  Section  [372]) .  By  fitting  a  broken 
powerlaw  function,  we  obtain  a  decay  index  of  —0.03  ±  0.03 
before  the  break,  which  steepens  to  —0.06  ±  0.02  thereafter. 
However,  further  observations  are  required  to  conhrm  that 
a  break  has  indeed  occurred. 

The  decay  parameters  that  we  find  for  EXO  0748-676 
are  comparable  to  that  obtained  bvlFridriksson  et  al.l  (  20101 
for  XTE  J1701-462.  These  authors  found  that  the  quies¬ 
cent  lightcurve  breaks  ~  20  —  150  days  post-outburst  and 
report  decay  indices  of  ~  —0.03  and  ~  —0.07  before  and 
after  the  break,  respectively.  iFridriksson  et  all  (20101  note 
that  possible  cross-calibration  effects  between  Chandra  and 
XMM-Newton  might  introduce  small  shifts  that  also  allow 
a  single  powerlaw  decay  with  slope  ~  —0.05.  The  cooling 
curves  of  KS  1731-260  and  MXB  1659-29  appear  to  have 
steeper  decays  with  indices  of  ~  —0.12  and  ~  —0.33,  re¬ 
spectively  (ICackett  et  al.|[2008il.  Due  to  the  scarcity  of  data 
points  it  is  unclear  whether  a  break  occurred  in  the  quies¬ 
cent  lightcurves  of  those  two  sources  (ICackett  et  al.ll2008l: 

I  Brown  fe  Cummindl2009lf . 

The  powerlaw  fits  show  no  indications  that  the  quies¬ 
cent  lightcurve  of  EXO  0748-676  is  levelling  off.  Thus,  it  is 
also  possible  that  the  neutron  star  temperature  continues  to 
decay  further  and  that  the  core  is  cooler  than  suggested  by 
the  exponential  decay  fits.  The  relatively  slow  decrease  of 
EXO  0748-676  might  then  reflect  that  the  crust  has  a  high 
conductivity,  albeit  lower  than  that  of  the  neutron  stars  in 
KS  1731-260  and  MXB  1659-29.  Further  observations  are 
thus  required  to  determine  whether  the  neutron  star  crust 
in  EXO  0748-676  has  nearly  cooled  down  and  to  be  able  to 
draw  conclusions  on  the  crust  and  core  properties. 
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Figure  6.  The  effective  temperatures  of  KS  1731-260  (green  di¬ 
amonds;  from  [Cackett  et  al  .11200(1.  MXB  1659-29  (red  bullets; 
fromlCackett  et  alJl200fl  l200fill.  XTE  J1701-462  (grey  crosses; 
from  IFridriksson  et  alj|20icf)  and  EXO  0748-676  (black).  Expo¬ 
nential  decay  fits  to  the  data  of  KS  1731-260,  MXB  1659—29  and 
XTE  J1701-462  are  shown  to  guide  the  eye  (green  dashed,  red 
dashed-dotted  and  grey  dotted  line,  respectively).  The  two  data 
points  of  XTE  ,11701-462  that  lie  above  the  decay  curve  are  likely 
due  to  a  temporal  increase  in  the  accretion  rate  onto  the  neutron 
star,  causing  reheating  (Fridriksson  et  al.  20 1  (I) . 
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